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Abstract. The origin of the far-infrared emission from the nearby radio galaxy M87 remains 
a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust 
emission, whereas others propose that the far-infrared emission can be explained by synchrotron 
emission without the need for an additional dust emission component. We observed M87 with 
PACS and SPIRE as part of the Herschel Virgo Cluster Survey (HeViCS). We compare the 
new Herschel data with a synchrotron model based on infrared, submm and radio data to 
investigate the origin of the far-infrared emission. We find that both the integrated SED and 
the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds 
with previous claims, we find no evidence of a diffuse dust component in M87. 
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1. Introduction 

At a distance of 16.7 Mpc, M87 is the dominant galaxy of the Virgo Cluster. It is 
one of the nearest radio galaxies and was the first extragalactic X-ray source to be 
identified. Because of its proximity, many interesting astrophysical phenomena can be 
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Figure 1. Left: the Herschel view of the central regions of M87. The bottom right image is a 
VLA 20 cm image from the FIRST survey. The 20 cm radio contours have been overlaid on the 
Herschel images. The field of view of all images is 160" x 90", beam sizes are indicated in the 
bottom right corner. Top right: the global SED of M87 from mid-infrared to radio wavelengths. 
Where no error bars are seen, they are smaller than the symbol size. The solid line in the plot is 
the best-fit power law of the ISOCAM, IRAS, MIPS, SCUBA, GBT, WMAP, and VLA data; the 
dotted line has only been fitted to the SCUBA, GBT, WMAP, and VLA data. Bottom right: 
residual between data and the best-fit synchrotron model in the infrared-submm wavelength 
range. The cyan line is a modified black-body model with T — 23 K and Ma = 7x 10 4 M©. 



studied in more detail in M87 than in other comparable objects. Particularly remarkable 
is the prominent jet extending from the nucleus, visible throughout the electromagnetic 
spectrum. The central regions of M87, in particular the structure of the jet, have been 



studied and compared intensively at radio, optical, and X-ray wavelengths, (e.g. Biretta 
et al.|1991MMeisenheimer et al || 19961 |Bohringer et al.|[200T| |Perlman et al.|2001| |Sparks 



et al.ll2004|lPerlman fc Wilson||2005| jKovalev et al.||2007| |Simionescu et al.j2008j ). 



Compared to the available information at these wavelengths, our knowledge of M87 at 
far-infrared (FIR) wavelengths is rather poor. A controversial issue is the origin of the 
FIR emission in M87, i.e., the question of whether the FIR emission is caused entirely by 
synchrotron emission or whether there is an additional contribution from dust associated 
with either the global interstellar medium or a nuclear dust component. Several papers, 
based on IRAS, ISO and Spitzer observations, arrive at different conclusion s (|Perlman| 
et al.||2007| |Buson et al]|2009| |Xilouris et aT|2004| |Shi et al.||2007| |Tan et al.|[2008| ). 



We investigate the nature of the FIR emission of M87 using new FIR data from the 
Herschel Space Observatory (Pilb ratt et al.|20f0 ), obtained as part of the science demon- 



stration phase (SDP) observations of the Herschel Virgo Cluster Survey (HeViCS,|Davies 



et al. 



2010). 



HeViCS is an approved open time key program, which has been awarded 
286 h of observing time in parallel mode with the PACS and SPIRE instruments. We 
will ultimately map four 4x4 square degree regions of the cluster at 100, 160, 250, 350 
and 500 /im, down to the 250 fim confusion limit of about 1 MJy sr _1 . The HeViCS 
SDP observations consisted of a single cross-scan of one 4x4 deg 2 field at the centre of 
the Virgo Cluster. While these SDP observations comprise only 6% of the total HeViCS 
observations, the analysis based on these observations already gives a prelude to the pri- 
mary HeViCS science goals including the effects of the environment on the dust medium 
of galaxies, the FIR luminosity function, the complete SEDs of galaxies and a detailed 
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Figure 2. A comparison between the synchrotron model image and the observed image at 
250 /im. The left panel shows the synchrotron image at the model resolution, the central panel 
shows the same model convolved to the SPIRE 250 /im beam and pixel size. The right panel 
shows the observed SPIRE 250 /im image. In all panels, the green lines are the contours of the 
synchrotron model at the model resolution. 



analysis of the dust content of dwarf and early- type galaxies ( Davies et al.|201Q Cortese 
eTaT|[20T0l |Clemens et aTJ|2010| |Smith et afl|Mo| |Grossi et al.||2010| |De Looze et al. 
2010| |Boselli et al.p OlO). The observations also enabled us to study the intensity and 
nature of the FIR emission of M87 flBaes et al.||2010). 



2. Analysis and conclusion 

The left part of Figure [l] shows the PACS and SPIRE images of the central 160" x 
90" region of M87, which is clearly detected in all five bands. The top right panel in 
Figure [I] shows the integrated SED in the infrared-submm-radio region between 15 /im 
and 100 cm, with the new Herschel data as well as ISOCAM, IRAS, MIPS, and SCUBA, 



GBT, WMAP, and VLA data gathered from the literature (|Xilouris et al. 2004 |Golombek 


et al.T988 Shi et al.|2007 Haas et al.|2004 


Cotton et al.|2009 Wright et al.|2009); the 


actual flux densities can be found in Baes et al. 


(2010). The solid line is the best-fit power 



law for all literature data and has a slope a = —0.76; the dotted line fits only the submm 
and radio data and has a slope a = —0.74. The bottom right panel in Figure [I] shows the 
residual from the best-fit power law in the infrared-submm wavelength region; clearly, 
the integrated Herschel fluxes are in full agreement with synchrotron radiation. The cyan 
line in this figure is a modified black-body fitted with T = 23 K and = 7 x 10 4 M@. 
This temperature is the mean dust equilibrium temperature in the interstellar radiation 
field of M87, determined using the SKIRT radiative transfer code (Ba es et al.|2003||2005[ ) 
and based on the photometry from Kormendy et al. (2009). The dust mass was adjusted 



to fit the upper limits of the residuals. It is clear that the SED of M87 is incompatible 
with dust masses higher than 10 5 M@. 

Although indicative, the analysis of the integrated SED does not definitively identify 
the origin of the FIR emission in M87. Approximating the global SED as a single power- 
law synchrotron model is indeed an oversimplification of the complicated structure of 
M87. Several studies have shown that M87 contains three distinct regions of significant 
synchrotron emission, each with their own spectral indices: the nucleus, the jet and asso- 
ciated lobes in the NW region, and the SE lobes (e.g., |Biretta et al.|1991 |Meisenheimer et 
al.|1996 Perlman et al.||2001 Shi et aL||2007 ). We have constructed a synchrotron model 
for the central regions of M87, based on a newly reduced MIPS 24 /im map and archival 



MUSTANG 90 GHz and 15, 8.2, 4.9, 1.6, and 0.3 GHz maps (Cotton et al. 2009). We 



fitted a second-order polynomial synchrotron model to each pixel of the MIPS + radio 
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data cube and used this synchrotron model to predict the emission of M87 at 250 /im 
(the SPIRE 250 jam image provides the optimal compromise between S/N and spatial 
resolution). Figure [2] shows the comparison between the synchrotron model prediction at 
250 /im and the SPIRE observations. At the model resolution, the three distinct compo- 
nents are visible, but when we convolve this synchrotron model image with the SPIRE 
250 jam beam, the three different components merge into a single extended structure with 
one elongated peak slightly west of the nucleus. Comparing the central and right panels 
of Figure |2j we see that the synchrotron model is capable of explaining the observed 
SPIRE 250 /im image satisfactorily. 

We conclude that for both the integrated SED and the SPIRE 250 jam map, we have 
found that synchrotron emission is an adequate explanation of the FIR emission. We do 
not detect a FIR excess that cannot be explained by the synchrotron model. In particular, 
we have no reason to invoke the presence of smooth dust emission associated with the 
galaxy interstellar medium, as advocated by Shi et al. (2007). For a dust temperature of 
23 K, which is the expected equilibrium temperature in the interstellar radiation field of 
M87, we find an upper limit to the dust mass of 7 x 10 4 M . Our conclusion is that, seen 
from the FIR point of view, M87 is a passive object with a central radio source emitting 
synchrotron emission, without a substantial diffuse dust component. 
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